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ABSTRACT 

The E/SO sequence of a cluster defines a boundary redward of which a reliable weak lensing signal can be 
obtained from background galaxies, uncontaminated by cluster members. For bluer colors, both background 
and cluster galaxies are present, reducing the average distortion signal by the proportion of unlensed cluster 
members. In deep Subaru and HST/ACS images of A1689, we show that the tangential distortion of galaxies 
with bluer colors falls rapidly toward the cluster center relative to the background reference level provided 
by the red background. We use this dilution effect to derive the cluster light profile and luminosity function 
to large radius, with the advantage that no subtraction of far-field background counts is required. The light 
profile of A 1689 is found to decline steadily to the limit of the data, r < 2 /i _1 Mpc, with a constant slope, 
d\og{L) / d\og{r) =-1.12±0.06, unlike the lensing mass profile which steepens continuously with radius, so 
that M/L peaks at an intermediate radius, ~ 100 h~ l kpc. A flatter behaviour is found for the more physically 
meaningful ratio of dark-matter to stellar-matter, when account is made of the color-mass relation of cluster 
members. We derive a cluster luminosity function with a flat faint-end slope of a = -1.05 ±0.07, nearly 
independent of radius and with no faint upturn to < -12. We establish that the very bluest objects are 
negligibly contaminated by the cluster ([V -i']ab < 0.2), because their distortion profile rises continuously 
towards the center following the red background ([V - i'~\ab > 1 -2), but is offset higher by ~ 20%. This larger 
amplitude is consistent with the greater estimated depth of the faint blue galaxies (z) ~ 2 compared to (z) ~ 0.85 
for the red background, a purely geometric effect. With a larger sample of background galaxies behind several 
clusters we may use this geometric effect to constrain the cosmological parameters in a model independent 
way, by comparing the weak lensing strength over a wide range of photometrically derived redshifts. Finally, 
we improve upon our earlier modeling of the mass profile by combining both the red and blue background 
populations, and very clearly exclude low concentration profiles predicted for massive CDM dominated haloes. 

Subject headings: cosmology: observations - gravitational lensing - galaxies: clusters: individual(Abell 1689) 
- galaxies: luminosity function 



1. INTRODUCTION 

The influence of "Dark matter" is strikingly evident in the 
centers of massive galaxy clusters, where large velocity dis- 
persions are measured and giant arcs are often visible. Cluster 
masses may be estimated by several means, leading to ex- 
ceptionally hig h central mass-to-lig ht ratios, M/Lr ~ 100 — 
300h(M/L R ) Q (ICarlberg et alJl200lh . far exceeding both the 
mass of stars comprising the light of the cluster galaxies and 
the mass of plasma derived from X-ray emission and the SZ 
effect. Reasonable consistency is claimed between dynam- 
ical, hydrodynamical and lensing-based estimates of cluster 
masses, supporting the conventional understanding of gravity. 
However, the high mass-to-light ratio requires that an uncon- 
ventional non-baryonic dark material dominates the mass of 
clusters, whose origin remains very unclear. 

In detail, discrepancies are often reported between masses 
derived from strong lensing and X-ray measurements, with 
the claimed X-ray masses often lower in the centers of clus- 
ters. High resolution X-ray emission and temperature maps 
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reveal that the majority of local clusters undergo repeated 
merging with s ub-clumps, and obvious shock fronts are see n 
in some cases dMarkevitch et alj|2002t iReiprich et al.ll2004t) . 
The double cluster IE 0657-56 (a.k.a., the "bullet" cluster, 
Z = 0.296) is the most extreme example studied, where the 
associated gas forms a flattened luminous shock-heated struc- 
ture lying between the two large distinct clusters, clearly in- 
dicating th ese two bodies collide d recently at a high relative 
velocity (Mark evitch et al.ll2004t) . with the gas remaining in 
between while the cluster galaxies have passed through each 
other relatively collisionlessly. Very interestingly, the weak 
lensing signal follows the double structure of the clusters, 
rather than the gas in between, cleanly demons trating that 
the bu lk of the matter is collisionless and dark dClowe et all 
2004). This system favors the standard cold dark matter 
(CDM) scenario, places a restrictive limit on the interaction 
cross-section of any fermionic dark matter, and disfavors a 
class of alternative gravity theories in which only baryons are 
present dClowe et al.l 120041) . Smaller but significant discrep- 
ancies o f the same sort are also claimed for o ther interacting 
clusters jNataraian et al.ll2002l:lJee et aT]|2005l) . 

Many clusters show no apparent signs of significant ongo- 
ing interaction; these have centrally s ymme t ric X- ray emis- 
sion and little obvious substructure dAllenl 119981) . and for 
some of these the lensing and X-ray (or dynamical) derived 
masses are claimed to agree, as would be expected for re- 
laxed systems (e.g., Arabadjis et al. 2002|; iRines et aT]|2003b 
iDiaferio et al.l 12005). More generally, since the dyanmics 
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of dark matter and most cluster galaxies is essentially colli- 
sionless, we would expect them to have similar radial pro- 
files. Biasing inherent in h ierarchical growth may signifi- 
cantly modify this similarity (Kauffman n et all 19971) . and dy- 
namical friction is expected to concentrate the relatively more 
massive galaxies in the core. These together with tidal effects 
may help to account for the unique properties of cD galaxies. 
Hence comparisons of the mass profile with the light profile 
of the cluster galaxies are expected to provide an additional 
insight into the formation of clusters and the nature of dark 
matter. 

Recent improvements in the quality of data useful for weak 
and strong lensing studies now allow the construction of much 
more definitive mass profiles that are sufficiently precise to 
test the predictions of popular models, relatively free of major 
assumptions. The inner mass profiles of several clusters have 
been constrained in some d etail via lensing, using multiply 
lensed background galaxies dKneib et aljl99 6: Hamm er" et al j 
19971: ISand et al Jl2002t jGavazzi et a l. 2003; Sand et al.ll20'04t 
Broadhurst et al.ll2005at ISharon et al.ll2005h . The statistical 
effects of weak lensing have been used to extend the mass 
profile to larger radii. The mass profiles derived from these 
observations have been claimed to show NFW-like behaviour 
dNavarro. Frenk. & White! I"l997lh with a continuous flatten- 
ing towards the center, but with higher density concentrations 
than expected. This is particularly evident for A 1689, for 
which we have constructed the highest quality lensing based 
mass profile to date, combining over 100 multiply-lensed im- 
ages and wea k lensing effects of disto rtion and magnification 
from Subaru (iBroadhurst et al .112005 allbh . 

A lesson learned from this earlier work is the importance 
of carefully selecting a background population to avoid con- 
tamination by the lensing cluster. It is not enough to sim- 
ply exclude a narrow band containing the obvious E/SO se- 
quence, following common practice, because the lensing sig- 
nal of the remainder is found to fall rapidly towards the clus- 
ter center, in contrast to the uncontaminated population of 
background galaxies lyi ng redward of the cluster sequence 
(IBroadhurst et al.ll2005bb . In the well studied case of A1689, 
there has been a long standing discrepancy between the strong 
and weak lensing effects, with the weak lensing signal under- 
predicting the observed Einstein radius by a fac tor of ~ 2.5 
(IClowe & Schneider! 1200 it iBardeau et afll2005h . based only 
on a minimal rejection of obvious cluster members using one 
or two-band photometry. 

In our recent weak lensing analysis of Subaru images of 
A 1689 the above behaviour was found when we rejected only 
the cluster sequence in the same way as others. This resulted 
in a relatively shallow trend of the weak lensing signal with 
radius and consequently an underprediction of the Einstein ra- 
dius ac counting for the earlier discrepancy (Broadhur st et al.l 
I2005bl) . If, however, one selects only objects redder than 
the cluster sequence for the lensing analysis, then the weak 
tangential distortion continues to rise all the way to the Ein- 
stein radius in very good agreement with the strong lensing 
strength. This red population is naturally expected to com- 
prise only background galaxies, made redder by relatively 
large k-corrections and with negligible contamination by clus- 
ter members, since the bulk of the reddest cluster members are 
the early-type galaxies defined by the cluster sequence. How- 
ever, for galaxies with colors bluer than the cluster sequence, 
cluster members will be present along with background galax- 
ies since the cluster population extends to bluer colors of the 
later-type members, overlapping in color with the blue back- 



ground. The effect of the cluster members is simply to reduce 
the strength of the weak lensing signal when averaged over 
a statistical sample, in proportion to the fraction of cluster 
members whose orientations are randomly distributed, there- 
fore diluting the lensing signal relative to the reference back- 
ground level derived from the red background population. 

We can turn this dilution effect to our advantage and use 
it to derive properties of the cluster population, in particular 
the radial light profile, for comparison with the dark matter 
profile. Deriving a light profile this way has advantages over 
the usual approach to defining cluster membership. The in- 
herent fluctuations in the number counts of the background 
population are a significant source of uncertainty in the usual 
approach of subtracting t he far-field level when defining the 
cluster population (e.g ., iPaolillo et alj 1200 U lAndreon et al.l 
120051: iPracy et aT1l2005l) . This uncertainty is often cited as a 
potential explanation for the substantial variation reported be- 
tween luminosity functions derived for different clusters, par- 
ticularly in the outskirts of clusters, where not only is the den- 
sity of galaxies lower, but their colors are bluer, thus harder to 
distinguish from the background using photometry alone. 

In §|2] we describe the observations and photometry of the 
Subaru images of A1689. In § [3] we describe the distortion 
analysis applied to the Subaru data. The distortion analysis of 
ACS images of A1689 is described in §2] In §|5]we describe 
the photometric redshift a nalysis of th e Suba ru and ACS im- 
ages with reference to the ICapak et al.l ([2004) sample of deep 
multi-color Subaru images. Our weak lensing dilution anal- 
ysis of the Subaru images is described in § [6] In § [7] we go 
on to derive the cluster luminosity profile and color, and in 
§[8]the cluster luminosity function is deduced at several radial 
positions. In § [9] we determine the M/L profile, and in § [10] 
we do a consistency check for the mass derived in this paper 
with previous estimations. Our conclusions are summarized 
in Slfll 

The concordance ACDM cosmology is adopted (O^ = 0.3, 
CI a = 0.7 but h is left in units of H Q /\QQ km s" 1 Mpc" 1 , for 
easier comparison with earlier work). 

2. SUBARU IMAGING REDUCTION AND SAMPLE SELECTION 

We have retrieved Suprime-Cam imaging of A 1689 in 

V (1920s) and SDSS i 1 (2640s) from the Suba ru archive 
SMOK A. 7 Reduction software developed by lYagi et aD 
d2002l) is used for flat-fielding, instrumental distortion correc- 
tion, differential refraction, PSF matching, sky subtraction, 
and stacking. The resulting FWHM is 0".82 in V and 0".88 
in i' with 0".202 pix" 1 , covering a field of 30' x 25'. 

Phot ometry is based on a com bined V + i' image using SEx- 
tractor dBertin & Arn outs 1996). The limiting magnitudes are 

V = 26.5 and i' = 25.9 for a 3a detection within a 2"aperture. 
We define three galaxy samples according to color and mag- 
nitude - "red", "green", and "blue" (see table[T]for summary), 
and for all our samples we define a limiting magnitude of 
i' < 26.5, to avoid incompleteness, as shown in FigureQ] The 
red galaxy sample consists of galaxies 0'".2 redder than the 
E/SO sequence of the cluster, which is accurately defined by 
the linear relation (V-i%/ S o = -0.03525/'+ 1.505, and up 
to 2 m .5 redder than this line to include the majority of the 
background red population. Very red dropout galaxies may 
be detected beyond this point. Indeed, one spectroscopically 
confirm ed example at g = 4.82 has been detected behind this 
cluster dFrve et al.ll2002l) . and such cases are excluded by this 

7 http://smoka.nao.ac.jp. 
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TABLE 1 
Sample Selection 



Sample Name 


mag 


limit 




ador limit 


N 


n [h 2 Mpc- 2 ] 


< z > 


<D> 


red 


18 < i' 


<26.5 


0.2 <(V- 


i')-(V-i') E/S o <2.5 


11774 


966.1 


0.871±0.045 


0.693±0.012 


green 


16 < V 


<26.5 


-0.3 < (V- 


-i')-(y-i') E/ so<0-t 


11963 


981.6 


1.429±0.093 


0.728±0.015 


blue 


23 < i' 


<26.5 


-1.2 < ev- 


')-(V-i') E/so <-0.45 


2459 


201.8 


2.012±0.124 


0.830±0.011 
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FIG. 1. — Color vs. magnitude diagram for A1689 cluster galaxies. The 
E/SO sequence is apparent at (V — i') ~ 0.8 where there is an overdensity of 
bright galaxies. The red points represent the background sample of galax- 
ies redder than the E/SO sequence. The blue points represent a background 
sample of the bluest objects in the field. The green points cover a range of 
color chosen to include the cluster sequence and bluer cluster members, but in 
addition background galaxies are also present whose colors fall in this range. 

upper limit, so that we do not need to make an uncertain cor- 
rection for the level of their weak lensing signal which will 
be significantly larger than for the bulk of the background red 
galaxy population. As we will show in § [5] most of these 
red background galaxies are at a much lower mean redshift 
of (z) ~ 0.85. The red galaxy sample is redder than the clus- 
ter sequence, made so by relatively large k-corrections, being 
largely comprised of early to mid-type galaxies at moderate 
redshift (see § [5j . Cluster members are not expected to extend 
to these colors in any significant numbers because the intrinsi- 
cally reddest class of cluster galaxies, E/SO galaxies, are de- 
fined by the cluster sequence and lie comfortably blueward of 
the chosen sample limit (see Fig.Q]), so that even large photo- 
metric errors will not carry them into the red sample. This can 
be demonstrated readily, as shown in Figure [2] where we plot 
the mean lensing strength as a function of color by moving the 
lower color limit progressively blueward, finding a sharp drop 
in the lensing signal at our limit, (V — i') < (V -i') E / so + 0.2, 
when the cluster sequence starts to contribute significantly, 
thereby reducing the mean lensing signal. 

We define the blue galaxy sample as objects 0"'.45 bluer 
than the sequence line, with the magnitude limit in the in- 
terval 23 < i' < 26.5, so as to take only the very faint blue 
galaxies which - as we establish below - are also negligibly 
contaminated by the cluster, with a weak lensing signal which 
has the same radial dependence as the red galaxy sample. We 
have explored the definition of the blue sample when we real- 
ized that the bluest objects in the field have a continuously ris- 
ing weak lensing signal towards the cluster center like the red 
galaxies, so that the "contamination" is minimal with an in- 
significant effect on the quantities of interest for our purposes. 
Figure|2]shows that as the blue sample upper color limit is ad- 
vanced redwards, the integrated strength of the mean weak 
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FIG. 2. — To establish the boundaries of the color distribution free of cluster 
members we calculate the mean tangential distortion averaged over the full 
radial extent of the cluster, done separately for the blue and red samples. On 
the right the red curve shows that gr drops rapidly when the bluer limit of 
the entire red sample is decreased below a color while lies +0"'.2 magnitudes 
redward of the cluster sequence. This sharp decline marks the point at which 
the red sample encroaches on the E/SO sequence of the cluster. On the blue 
side the red limit of the blue sample is chosen to lie —0'" .45 blueward of the 
cluster sequence, marking the point at which significant contamination by 
the cluster acts to dilute the weak lensing signal from later-type bluer cluster 
members. 

lensing signal declines markedly within 0'".45 of the cluster 
sequence, at (V — i ! ) > (V-/')£/so - 0.45. The reduction in 
signal is more gradual than for the red population (both illus- 
trated in Fig. |3 because the blue cluster members do not lie 
along a sharp sequence but contribute a diminishing fraction 
relative to the background at bluer colors. 

The green galaxy sample is simply selected to lie be- 
tween the red and blue samples defined above, (V - j')e /so*q I 
(Fig. H}, with generous limits set to include the vast major- 
ity of cluster galaxies, since - as we have established - both 
the red and blue samples are negligibly contaminated by clus- 
ter members, and hence the vast majority of cluster members 
must lie within this intermediate range of color. A narrow gap 
on each side of these samples is left out of our analysis to en- 
sure that the definition of the background does not encrouch 
on the cluster population. Note that unlike the green sample 
containing the cluster population, the background populations 
do not need to be complete in any sense but should simply 
be well defined and contain only background. Increasing the 
green sample to cover these narrow gaps does not lead to any 
particularly significant change in our conclusions, but only 
increases somewhat the level of noise by including relatively 
more background galaxies. Within the green sample there are 
of course background galaxies, and the purpose of this pa- 
per is to make use of the relative proportion of these cluster 
and background populations via weak lensing to establish the 
properties of the cluster galaxy population, by using the the 
dilution of the weak lensing signal of the background galax- 
ies due to the cluster members. 
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3. DISTORTION ANALYSIS OF SUBARU IMAGES 

We use the IMC AT package developed by N. Kaiser 
8 to perform object detection, photometry and shape 
measurements, following th e form alism outlined in 
dKaiser. Squires. & Broadhurstl I1995L hereafter KSB). 
We have modified th e method somewh at following the 
procedures described in Erb en et all d200U see Section 5). 

To obtain an estimate of the reduced shear, g a = 7 a /(l — 
k), we measure the image ellipticity e a from the weighted 
quadrupole moments of the surface brightness of individual 
galaxies. Firstly the PSF anisotropy needs to be corrected us- 
ing the star images as references: 

e' a =e a -P«% (1) 
where P sm is the smear polarizability tensor being close to 
diagonal, and q* a = (P* sm )^gS* * s tr, e stellar anisotropy ker- 
nel. We select bright, unsaturated foreground stars identi- 
fied in a branch of the half-light radius (r/,) vs. magnitude 
(i 1 ) diagram (20 < i' < 22.5, (r A ) median = 2.38 pixels) to cal- 
culate q* a . In order to obtain a smooth map of q* a which is 
used in equation (HJ, we divided the 9K x 7.4K image into 
5x4 chunks each with 1.8K x 1.85K pixels, and then fit- 
ted the q* in each chunk independently with second-order bi- 
polynomials, q"(0), in conjunction with iterative rr-clipping 
rejection on each component of the residual e* — P*Snq*a{&). 
The final stellar sample consists of 540 stars, or the mean sur- 
face number density of n* = 0.72 arcmin -2 . From the rest of 
the object catalog, we select objects with 2.4 < r/, < 15 pixels 
as an i'-selected weak lensing galaxy sample, which contains 
61, 115 galaxies or n g ~ 81 arcmin" 2 . It is worth noting that 
the mean stellar ellipticity before correction is (e~i*,e~2*) — 
(-0.013,-0.018) over the data field, while the residual e* 
after correction is reduced to e™ s = (0.47 ± 1.32) x 10" 4 , 
e* res = (0.54 ± 0.94) x 10" 4 . The mean offset from the null ex- 
pectation is |e* res | = (0.71 ±1.12) x 10" 4 . On the other hand, 
the rms value of stellar ellipticities, a„ = (|e*| 2 ), is reduced 
from 2.64% to 0.38% when applying the anisotropic PSF cor- 
rection. 

Second, we need to correct the isotropic smearing effect on 
image ellipticities caused by seeing and the window function 
used for the shape measurements. The pre-seeing reduced 
shear g a can be estimated from 

g a = (Pg l )*pe'i3 (2) 

with the pre-seeing shear polarizability tensor P£g. We fol- 
low the procedure de scribed inlErben et al. (2001) to measure 
P 8 (see also § 3.4 of iHetterscheidt et al.ll2006l) . We adopt the 
scalar correction scheme, namely, 

P s a0 =\^[P s ]S a = P^ a (3 (3) 

dHudson et al.lll998t IHoekstra et al.|[l99l lErben et~aT]|2001b 
IHetterscheidt et al.l 12006b . The P? measured for individual 
objects are still noisy especially for small and faint objects. 
We thus ado pt a smoothing scheme in ob ject para meter space 
proposed bvlVan Waerbeke et al.l d2000l see also lErben et ail 
l200UlHamanaet al.ll2003l) . We first identify thirty neighbors 
for each object in r g -i' parameter space. We then calculate 
over the local ensemble the median value (P|) of P* and the 
variance er 2 of g = g\ + ig2 using equation ©. The dispersion 

8 http://www.ifa.hawaii/kaiser/IMCAT 



a g is used as an rms error of the shear estimate for individual 
galaxies. The mean variance ct 2 over the sample is obtained 

as ~ 0.171, or y^ 2 - 0.41. 

In the previous study by Broadhurst et al. (2005b), those 
objects that yield a negative value of the raw PJ estimate were 
removed from the final galaxy catalog to avoid noisy shear 
estimates. On the other hand, in the present study, we use all 
of the galaxies in our weak lensing sample including galax- 
ies with P' < 0. After smoothing P£ in the object parameter 
space, all of the objects yield positive values of (Pf), with the 

minimum of as 0.04. The median value of (Pi) over the weak 
lensing galaxy sample, including galaxies with P! < 0, is cal- 
culated as w 0.32. For a reference, the sub-sample of galax- 
ies with P! > gives the median average of w 0.33, mostly 
weighted by galaxies with r g = 2-2.5 pixels. Finally, we use 
the following estimator for the reduced shear: 

g a =e'J(Pl)- (4) 

The quadratic shape distortion of an object is described by 
the complex reduced-shear, g = g\+ igi- The tangential com- 
ponent gj is used to obtain the azimuthally averaged distor- 
tion due to lensing, and computed from the distortion coeffi- 
cients gi,g 2 : 

g T = -(g\cos26 + g2sm20), (5) 

where 8 is the position angle of an object with respect to the 
cluster center, and the uncertainty in the gj measurement is 
a = a g /v2 in terms of the rms error a g for the complex 
shear measurement. The cluster center is well determined 
from sy mmetry of the strong lensing pattern (Broadhur st et al.l 
2005a). The estimation of gj only has significance when eval- 
uated statistically over large number of galaxies, since galax- 
ies themselves are not round objects but have a wide spread in 
intrinsic shapes and orientations. In radial bins we calculate 
the weighted average of the gjs and the weighted error: 

(g T (r„)) = ^ll4 ( 6 ) 

L 1 / cr 

VT(r„)=(j2 l /< j2 )~ 1/2 ■ (7) 

It has been shown that such weights depend on the size 
of the objects but m ostly on their magnitudes (see e.g. 
IHoekstra et alj d2006l) ). Therefore, as apparent magnitude 
increases with redshift the redshift distribution of sources 
will be modified to some extent by this weighti ng scheme. 
We ha ve investigated this using the catalogue of Cap ak et ail 
(2004) as here photometric redshifts are estimated (see § [5] 
for a fuller desrciption of the photometric properties of this 
sample). First, we generated from the Capak catalog blue/red 
background galaxy samples with the same color-magnitude 
criteria as the present study. We then derived an i' -magnitude 
vs. photo-z relation for each galaxy sub-sample. We subdi- 
vided the data into magnitude (i') bins and derived a magni- 
tude (i') vs. photo-z relation using median averaging. 

We then assume this magnitude-redshift relation holds in 
our A 1689 data and obtain for each galaxy in A 1689 an es- 
timate of redshift via the magnitude - photo-z relation. It is 
then straightforward to have an effective redshift distribution 
taking into account weak lensing statistical weights, w. We 
can see a qualitative feature that although low-z background 
galaxies are more strongly weighted than higher-z ones, the 
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FIG. 3. — Tangential shear profiles, gj(r) of the red and blue background 
populations. The tangential distortion profile of both decline smoothly from 
the center, remaining positive to the limit of the data. The red galaxies are 
fitted well with a simple power-law, d\oggj /d\og r = —1.17 ±0.1, The blue 
sample is more noisy but also well represented by the same relation, only 
offset in amplitude by 23 ± 17% and is related to the greater depth of the blue 
population relative to the red, §|3] 

effect on our observed redshift distribution is negligible be- 
cause our redshift selection window does not sample these 
larger angle, lower redshift objects, but rather the more dis- 
tant faint population whose small anglular sizes are heavily 
influenced by the seeing. 

In Figure |3] we compare the radial profile of gj of the red 
and blue samples defined above. These have a very similar 
form implying that the blue sample, like the red sample, is 
dominated by background galaxies with negligible dilution 
by cluster members even at small radius where the cluster 
overdensity is large. Very interestingly a clear offset is vis- 
ible between these profiles over the full range of radius, with 
the amplitude of the blue sample lying systematically above 
that of the red sample, as shown in Figure [3] This is read- 
ily explained as a depth related effect, as we show below in 
§|5] where we evaluate the redshift distributions of these two 
populations. We find the blue sample to be deeper than the 
red, and since lensing scales with depth, so should the lensing 
profiles be offset from each other by the same scale. 
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FIG. 4. — Tangential shear vs. radius. The green points represent the 
intermediate color sample, containing both cluster and background galax- 
ies. The black points show the level of tangential distortion of the combined 
red+blue sample of the background. The green points fall close to the back- 
ground level defined at large radii, indicating the green sample is dominated 
by background galaxies, and falls short towards the cluster center where clus- 
ter members increasingly dilute the lensing signal. 
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Fig. 5. — Color vs. magnitude diagram for the central region covered 
by ACS photometry with magnitudes transformed to AB system to match 
the V,I photometry from Subaru. Notice the prominent E/SO sequence and 
the greater depth of these data compared with Subaru shown in Figure 1 . A 
one-to -one comparison of magnitudes for objects in both datasets is shown in 
Figure fl5] 

The tangential distortion of the green population behaves 
quite differently (Fig. falling well below the background 
level near the cluster center. The green sample has a maxi- 
mum signal at intermediate radius, 3' -5', and then declines 
quickly inside this radius as the unlensed cluster galaxies 
dominate over the background in the center. Notice that the 
green sample does not fall to zero at the outskirts but rises up 
to almost meet the level of the background sample, indicating 
that the majority of the green sample at large radius comprises 
background rather than cluster members. We go on to use the 
ratio of the distortion of the green sample compared with the 
background level to determine the proportion of cluster mem- 
bers in § [6j but to do so we first evaluate the expected depths 
of our samples in § [5] in order to make a precise comparison 
of the lensing signals between them, since the lensing signal 
scales geometrically with increasing source distance and must 
be accounted for in any comparisons. 

The results of this paper depend on the ratio of the 
background distortion to the cluster contaminated distortion, 
{g~ T jgj ), so that the 5-10% level calibration correction fac- 
tors estimated from simulations done by the STEP project 
dHevmans et al.ll2006l) for the various weak distortion meth- 
ods are not of major concern for the bulk of our work. 

4. DISTORTION ANALYSIS OF ACS/HST IMAGES 

In the center of the cluster inside a radius of approximately 
1' the Subaru data become limited in depth by the extended 
bright haloes of the many luminous central galaxies. This 
region is far better resolved and more deeply imaged with 
HST/ACS in 20 orbits of imaging shared between the g'r'i'z' 
passbands. Many multiple images are known here, defining 
accurat ely the shape of both th e tangential and radial critical 
curves (Broadhurst et al. 2005a). Here we analyze the statis- 
tical distortion of the shapes of the many galaxies recorded 
in these images, to extend our analysis of the properties of 
the cluster galaxies into the center, allowing an accurately de- 
fined central cluster luminosity function to faint luminosities. 
In addition, it will be interesting to see how consistent the dis- 
tortion profile derived here independently matches the mass 
profile obtained previously from the strongly lensed multiple 
images. 

We stick to very similar definitions of the three-color se- 
lected population as with Subaru, but extend their depths by 
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an additional l m .5 magnitudes since the ACS data are so 
much deeper. The ACS images are limited to m=28.5 {5a) 
in each of the passbands. The reduction of the ACS image 
and th e photometry for the fai nt sources is described in de- 
tail in iBroadhurst et a l. (2005a), including the subtraction of 
the bright central galaxies in the cluster which is essential 
for obtaining accurate photometry and shape measurements 
of central lensing images including radial arcs and demagni- 
fied central images. For the dist ortion analysis we p refer the 
IM2SHAPE method developed bv lBridle et all d2002l) for deal- 
ing in particular with relatively elongated images produced by 
lensing in the strongly lensed region. This is an improvement 
over the standard KSB method which we used in the weak 
lensing regime appropriate for everything except the central 
region r < 2', and used for the Subaru analysis described 
above. 

Using this method, galaxies are fit to a sum of two sheared 
Gaussians convolved with a PSF. Each Gaussian has two free 
parameters, amplitude and width. The centroid and the shear 
are also allowed to vary, but these are restricted to be the same 
for both Gaussians. Meanwhile, the PSF fo r each g a laxy i s 
determined based on models described in iJee et al.l d2005l) . 
These PSF models for ACS's WFC were derived from ob- 
servations of the globular cluster 47 Tuc (PROP 9656, P.I. 
De Marchi). As the distortion measurements are performed 
in the detection image, each galaxy is assigned an "average" 
PSF based on the different filters and chip positions in which 
it was observed. 

We plot the resulting values of gr(r) (Fig. [6) for the blue, 
green and red galaxies defined in the same color ranges as the 
Subaru data, but to fainter magnitudes. A very well defined 
saw-tooth pattern is visible, showing that images are maxi- 
mally radially aligned at about 17" and then maximally tan- 
gentially aligned at about 47". This is a very clear signature 
of strong lensing, where the maximum corresponds to the lo- 
cation of the tangential critical curve (Einstein radius), and 
the minimum to the radial critical curve, where images are 
maximally stretched in the radial direction generating a ring 
of long images pointing to the center of mass, as found in 
IBroadh urst et al. (2005a). The location of these critical radii 
agrees very well with those derived from t he model to the 
strong lensing data for this cluster, fitted by Broadhur st et al.l 
d2005ah . 

Another clearly defined radius can also be identified from 
the point where the images distortion goes through zero, gj = 
0, at a radius in between these two critical radii at about r ~ 
27". It is important to note that in this region, between the 
two critical curves, the parity is p = — 1 (odd parity), and here 



1 



(8) 



(see iKaiserll 1995b . Here, instead of measuring g we are mea- 
suring — 1 < e < 1. (This is different than in the weak lensing 
region, outside the tangential critical curve, where p = +1 and 
g = e, and therefore no distinction needed to be made). Since 



' = — oc 1 — K, 



(9) 



zero distortion corresponds to re = 1 curve, which lies in be- 
tween the tangential and critical curves. 

Note, for several reasons we cannot expect that the data 
will reach the theoretically extreme value of gj = 1 at the tan- 
gential critical radius, meaning that the images are infinitely 
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FIG. 6. — Distortion profile for the central r < 2' area covered by ACS 
(filled circles) together with the weaker distortions measured by Subaru at 
larger radius (open symbols). The black points represent the combined blue 
and red sample of background galaxies, and the green points include the clus- 
ter members. A remarkable saw-tooth pattern is visible for the background 
galaxy distortions in the strong lensing region where the tangential and radial 
critical radii are clearly visible corresponding to a maximum and a minimum 
in the value of gj, respectively. In between the distortion passes through zero, 
where the degree of tangential and radial distortion is equal, leaving images 
unchanged in shape at a radius where n = 1 . The distortion of the green 
sample is consistent with zero in the inner region where the cluster members 
dominate the sample but at larger radius the green and black points merge 
for r > 3', indicating that there the green sample comprises predominately 
background galaxies. 

stretched tangentially, and also gj=-\ at the radial critical ra- 
dius where they are stretched infinitely in the radial direction. 
By definition, weak lensing measurements will underestimate 
the sting distortions near critical curves. In addition, convo- 
lution by the redshift distribution of the background sources 
will smooth these features out. Nonetheless, we can define 
their locations in radius rather precisely and these positions 
must be reproduced in any satisfactory model, at r ~ 17" and 
r ~ 48". In addition the radius at which g T = is also well 
defined at about r ~ 28" and corresponds to a surface density 
where k = 1 , supplying another important constraint on model 
mass profiles. 

Outside the tangential critical curve we find that the tangen- 
tial shear of the background sample (black circles in Fig. |6]l 
drops to gr.B ~ 0.2 at r = 2', in good agreement with the Sub- 
aru analysis at this radius, giving us confidence in the consis- 
tency of our work. For the color range defined above, which 
includes the cluster sequence and all the bluer members of the 
cluster galaxy population, gr.c ~ over the full range of ra- 
dius of the ACS data (green circles in Fig. [6}, indicating - as 
expected - that the galaxy population in the this color range 
is dominated by cluster members with negligible background 
contamination. 

5. PHOTOMETRIC REDSHIFTS 

We need to estimate the respective depths of our color- 
magnitude selected samples when estimating the cluster mass 
profile, because the lensing signal increases with source dis- 
tance, and therefore must differ between the samples. The 
effect of this difference in distance on the weak lensing signal 
is simply linear as we can see from the relation between the 
dimensionless surface mass density, 



where 



re(r) = S(r)/S crft , 



4-nGdi di s 



(10) 
(11) 



7 



3.5 - 
3 - 

2.5 - 



3 - 

1.5 - 




-1 1 ' ' 1 1 1 ' 1 

1B 18 20 22 24 26 26 

'«B 

FIG. 7. — Color-magnitude diagram for Capak galaxy catalog for the same 
passbands as the photometry of A 1689 shown in FigfT] This photometry is of 
a field region and serves as our reference for evaluating the expected depth of 
the background samples defined in §|2] 

and the tangential distortion: 

(gT(r)) = (R(r)-K(r))/(l-K(r)) (12) 
so that in the weak limit where k is small, 

(g r (r))cx^(S(r)-£(r)) (13) 

and hence for an individual cluster, with a fixed redshift and 
a given mass profile, the observed level of the weak distortion 
simply scales with the lensing distance ratio. Further details 
are presented in the appendix. The mean ratio of d\ s j d s , which 
is weighted by the redshift distribution of the background pop- 
ulation corresponding to our magnitude and color cuts, is cal- 
culated using the expression 

a I ^(z)N(z)dz 

Since we cannot derive complete samples of reliable pho- 
tometric redshifts from our limited 2-color V,i' images of 
A 1689, we instead make use of other deep field photometry 
covering a wider range of passbands, sufficient for photomet- 
ric redshift estimation of the faint field redshift distribution 
appropriate for samples with the same color and magnitude 
limits as our red, gree n and blue populat ions. 

The photometry of ICapak et"all (120041) is very well suited 
for our purposes, consisting of relatively deep multi-color 
photometry over a wide field taken with Subaru, producing re- 
liable photometric redshifts for th e majority of field g alaxies 
to faint limiting magnitudes. The ICapak et alj d2004l) galaxy 
catalog contains almost 50,000 galaxies over 0.2 sq. deg. 
with UBVRIZ photometry. We have estimated photometric 
redshift s for th is catalog using the Bayesian based method of 
iBenftezI ([2000), with a prior based on the redshift and spec- 
tral type distributions of th e HDF-N, with a spe ctral library 
containing the templates of iBenftez et alj d2004l) with an ad- 
ditiona l two blue starburst galaxies as described in lCoe et"aT] 
(120061) . 

A full redshift probability distribution is produced for each 
galaxy of the form: 

p(z\Q(xJ2p(z,T\m )p(C\z,T), (15) 

T 

where p{C\z, T) is the redshift likelihood obtained by com- 
paring the observed colors C with the redshifted library of 



templates T. The factor p(z,T\mo) is a prior which repre- 
sents the redshift/spectral mix distribution as a function of 
the observed /-band magnitude. We use a prior which de- 
scribes the redshift/spectral type mix in the HDF-N, which 
has been shown to significantly reduce the number of "catas- 
trop hic" errors (Az > 1) in the photometric redshift catalog 
(see IBenftez et al.l 12004 and references therein). For each 
galaxy we look at its redshift probability distribution p(z) and 
identify up to 3 redshift local maxima. Each of these maxima 
corresponds to a redshift Zi, spectral type t„ and a discretized 
probability p,(z,,f,) < 1. Using this information we generate 
a mock observation of all the zt,U combinations in the Subaru 
filters, and then build a redshift histogram by selecting galax- 
ies using the same color cuts and adding up their probabilities 
in each redshift bin. 

The color-magnitude diagram for the Capak catalog galax- 
ies is shown in Figure |7l where the equivalent color- 
magnitude selected samples are displayed. The resulting 
mean redshift of the background galaxies in each of our three 
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FIG. 8. — The weighted mean lensing depth di s /d s as a function of the 
apparent i' -magnitude limit of the red, green & blue backgrounds, calculated 
using the photometric redshifts of the Capak et al. 1 2004) sample. The ex- 
pected depth of the samples differs significantly between the samples and in 
general the distance ratio grows only slowly with increasing apparent magni- 
tude for each sample. 



0.85 - 

0.8 - 




0.7S r 




24.5 25 25.5 26 26.5 27 

Limiting magnitude 



FIG. 9. — Similar to the previous figure but the redshift is plotted against 
apparent magnitude corresponding to the mean depth calculated above. The 
average redshift differs significantly between the three samples, being lowest 
for the red background (z.) ~ 0.85 and highest for the blue (z) ~ 2. 
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color-selected samples is caclulated as a function of limiting 
magnuite of the sample (Fig. [9]), by using the redshift dis- 
tribution from the Capak catalog. The redshift distribution is 
also used to evaluate the weighted mean depths (D) (shown in 
Fig. [8] as a function of sample limiting magnitude), for com- 
paring the weak lensing amplitudes between the green and 
the red+blue samples. This is done by dividing up each sam- 
ple into 81 independent bins of 2', calculating the weighted 
mean redshift and depth in each, and taking the mean value 
and variance over the bins. The mean redshift of the red 
sample is only (z r ed) ~ 0.871 ± 0.045, whereas the blue sam- 
ple is calculated to have, (zbiue) ~ 2.012 ±0.124. The green 
sample lies in between with, (z gr een) ~ 1.429 ±0.093. The 
weighted relative depths of these samples using equation ( TBI , 
for samples selected to our magnitude limit of i < 26.5, are 
(D red ) = 0.693 ±0.012, (D gree „) =0.728±0.012, and (D blue ) = 
0.830 ±0.011, and the corresponding redshifts zd equiva- 
lent to these mean depths are, ZD,red = 0.68, Zo,green = 0.79, 
ZD,Mue = 1 -53, respectively. Hence the ratio of the mean depth 
of the blue sample to the red sample is (Di,i ue ) / (D re d) = 1 -20, 
accounting well for the observed offset seen in Figure[3] 

We also make use of the Capak "green" sample to inves- 
tigate the level of "cosmic variance" in di s /d s , and although 
there is variation in the redshift distribution the variance of 
the mean redshift is remarkably tight, and as quoted above 
we find a very small variance associated with the mean lens- 
ing depth, a{(di s /d s )) = 0.015. This stability is also a feature 
noticed in pencil beam redshift surveys in general, that the 
mean depth is stable to spikes in the redshift distribution, e.g., 
IBroadhurst etaTI (Il988l) . 

The form of the distance ratio D can be expressed in terms 
of the redshifts of the source and lens for a given set of cos- 
mological parameters. In the main case of interest, that of a 
flat model with a nonzero cosmological constant, the relation 
is given by 



d, s _ (l+z/) 2 Cfe) 

d s (l + Zs) 2 ((Zs) 



C(*) = 



dz 



[fi A + ft M (l+z)3]l/2- 



(16) 
(17) 



General expressions for the dependence of this distance ratio 
on arbitrary combinations of £7m a nd £7a are lengthy and can 
be found in Fukugit aet al.1 (Il990h . For a low redshift clus- 
ter like A1689 (z=0.183), t he form of this function is rather 
flat for sources at z > 1, see IBroadhurst et al.l d2005al) . There- 
fore the main uncertainty in determining the cosmological pa- 
rameters from a comparison of gr between samples of differ- 
ent redshifts, is small compared with clustering noise along a 
given line of sight behind the cluster, as examined in detail by 
(Broadhurs t et al.ir2005al) . Thus, we do not seriously examine 
this effect here but rather simply a dopt the recent (thre e year) 
WMAP cosmological parameters ( Spe rgel et al.l l2006) when 
making the above depth correction. With sufficient number of 
clusters and similar or better photometric redshift information 
(from multiple filter observations) one can hope to examine 
the trend of redshift vs. lensing distance in the future. 

Note that lensing magnification, \i, will modify gj slightly 
by increasing the depth to a fixed magnitude. But the mag- 
nification is small, /i < 0. m 2 over most of the cluster, r > 3'. 
In any case, the dependence of the mean redshift on depth 
is a slow function of redshift, so that it is safe for our main 
purposes to ignore the effect of magnification on the depth 
of our samples. Furthermore, since we are only interested in 
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FIG. 10. — Tangential distortion of the bright cluster sequence (i' < 21"'. 5) 
galaxies plotted against radius from the cluster center. By choosing the bright 
part of the sequence we minimize the background contamination, and can 
therefore check that the tangential distortion of the cluster members is negli- 
gible, which indeed is very clear from this figure. 

the proportion of the cluster relative to the background for 
our purposes, we are not affected by the modification of the 
background number counts caused by lensing, which has been 
shown to significantly deplete the surface density of back- 
ground red galaxies in A 1689, and found to be consistent with 
the predicted l evel of magnification ba sed on the distortion 
measurements (Broadhurst et al. 2005b). 

6. WEAK LENSING DILUTION 

We can now estimate the number density of cluster galax- 
ies by taking the ratio of the weak lensing signal between 
the green sample and the background sample, with the back- 
ground including both red and blue galaxies selected from the 
Subaru catalog as explained in §12 and accounting for differ- 
ences in the relative depths of these samples, as explained in 

§13 

Cluster members are unlensed and hence assumed to have 
random orientations, so they are expected to contribute no net 
tangential lensing signal. This assumption can be examined 
for the brighter (/' < 21'". 5) cluster sequence galaxies whose 
tight sequence protrudes beyond the faint field background 
(Fig. [B with negligible background contamination, so that 
we are secure in selecting this subset to test the assumption 
that the cluster galaxies are randomly oriented. Indeed, the 
net tangential signal of this population is consistent with zero 

(Fig. [[0]), gf = 0.0043 ±0.0091. 

For a given radial bin (r„) containing objects in the green 
sample, whose width in color has been chosen to encompass 
the full range of cluster galaxies (§0), the mean value of 
(eq. |6]) is an average over background and cluster members. 
Thus, its mean value (g^ ) will be lower than the true back- 
ground level denoted by (gj< ) (Fig. H} in proportion to the 
fraction of unlensed galaxies in the bin that lie in the cluster 
(rather than in the background), since the cluster members on 
average will add no net tangential signal. Therefore, 



fd(r n ) 



Ncrea 



L- 



(18) 



is the cluster membership fraction of the green sample (see 
full derivation in the appendix). 

Thus, we can use this effect to quantify statistically the 
number of cluster galaxies by comparing gj, ■* with the true 
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FIG. 1 1 . — Fraction of cluster membership vs. radius. Cluster membership 
is proportional to the dilution of the distortion signal of the green sample, 
relative to the expected distortion of the background galaxies set by the red 
and blue samples. As expected, close to the cluster center the fraction of 
cluster members in the green sample becomes maximal, whereas at larger 
radius r > 3', the fraction of cluster members is small, indicating that most 
of the green sample comprises background galaxies. 



background level derived from the pure background red and 
blue samples g ( j\ at a fixed radius. This is shown in Fig- 
ure[]Tl where we have also taken into account the effect of the 
relative depths of the differing samples. We find that the frac- 
tion of cluster members drops smoothly from ~ 100% within 
r < 2' to only — 20% at the limit of the data, r ~ 15'. 

7. CLUSTER LIGHT AND COLOR PROFILES 

To determine the luminosity profile of the cluster galax- 
ies, we need to go further, because in general the brightness 
distribution of the cluster members is different than that of 
the background galaxies; specifically, it is skewed to brighter 
magnitudes, certainly for the bulk of the cluster sequence. To 
account generally for any difference in the brightness distri- 
butions we can subtract a "^-weighted" luminosity contri- 
bution of each galaxy, which when averaged over the distri- 
bution will have zero contribution from the unlensed cluster 
members. We first calculate the "g-weighted" correction in 
arbitrary flux units. We estimate the total flux of the cluster in 
the n' h radial bin, 



F d (r„) = y>; 



(G) 



<Sr(r„)®> 



(19) 



where the sum is over all galaxies in the radial bin. 

The flux is then translated back to apparent magnitude, and 
from that the luminosity is derived. First we calculate the 
absolute magnitude, 



Mj, =i'-5\ogd L -K(z) + 5, 



(20) 



where the -correction is evaluated for each radial bin accord- 
ing to its V - i' color, which - after the correction is made for 
each of the bands - is now the cluster color. The luminosity is 
then 

L ., = 1O O.4(M, -M,) i;7Q; (21) 

where M,>q = 4.54 is the absolute i' magnitude of the Sun (AB 
system). 

The result yields the luminosity profile of the cluster as 
shown in Figure Q~2] (red squares). Here we can see that 
the cluster luminosity profile is well approximated by a sim- 
ple power-law with a projected slope of dlog(Lj>) / dlog(r) ~ 
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FIG. 12. — The "g-Weighted" luminosity density vs. cluster radius. Each 
galaxy's flux F{ (green sample) is weighted by its tangential distortion gj , 
with respect to the background distortion signal (red points). The filled circles 
represent the ACS data, and the empty squares are for the Subaru data. The 
blue points are derived from integrating over the luminosity functions of the 
same radial bins (see §[§), and serve as a consistency check, showing good 
agreement between these differing calculations. The dashed line is the best 
fitting linear relation for the blue points. 

-1.12 ±0.06, to the limit of the data. We also show in Fig- 
ure [13] the unweighted luminosity profile with no correction 
for the field, demonstrating that the g-weighted correction is 
negligible at small radius as expected, since the cluster domi- 
nates numerically over the background, but becomes increas- 
ingly more important at larger radius where the background 
dominates. Note that we derive a more accurate inner lumi- 
nosity profile using the ACS photometry for the central region 
(Fig- El re d circles), and here there is only a negligible cor- 
rection for the background due to the high central density of 
galaxies in this cluster. 

In a carefu l study of clu s ters a nd groups identified in the 
SDSS survey. [Hansen et all d2005l) find a similar slope for the 
most massive clusters, in terms of the composite surface den- 
sity profile of c/log(«)/iilog(r) ~ -1.05 ±0.04, over the ra- 
dius range r < 2Mpc, with slightly shallower slopes occur- 
ring in the less overdense clusters and groups. This may be 
compared directly with our slope derived above, assuming a 
constant M*/L, for the ratio of galaxy mass to galaxy lumi- 
nosity. More directly we derive a density profile using the 
/ C ;« (G) which also gives a slope of ~ -0.9 ± 0.09. 

In the same manner, we construct a "g-weighted" color pro- 
file of the cluster, V—i', which we have corrected as described 
above. We obtain a color profile that shows a weak tendency 
towards bluer colors with increasing radius, as expected, in- 
dicating a tendency towards later-type galaxies at large radius 
(Fig.[T4"ii. Also shown is the unweighted color profile (green 
points), which again is steeper due to the uncorrected field 
component which dominates numerically over the cluster at 
larger radius and is generally bluer in color than the cluster. 
This change in color with radius corresponds to a significant 
radial gradient in spectral type, from predominantly early- 
type with (V - i')ab — 0.84, to mid sequence type, Sb, with 
(V-i')AB — 0.73, and indicates that for this cluster very blue 
starburst and Scd galaxies are not the dominant population at 
the limiting radius of our sample (r ~ 2 /r'Mpc), where oth- 
erwise the col or would tend to (V- i')AB — 0.5, using standard 
template sets (iBenftez et a l. 2004). We go on to make use of 
this color-radius relation in § [9] when examining the radial 
profile of the ratio of total cluster mass to the stellar mass 
in galaxies. We do this by correcting the luminosity profile 
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FIG. 13. — The "g-Weighted" luminosity density vs. cluster radius (red 
squares), compared to the unweighted luminosity density (green circles), 
showing as expected the increasing size of the correction with increasing 
radius, where the sample becomes increasingly dominated by background 
galaxies. 
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FIG. 14. — Galaxy color profile after weighting the color of each object 
by its individual distortion, g;, accounting for any difference between the 
color distribution of the cluster and background populations comprising the 
green galaxy population. The color of the cluster members becomes slowly 
bluer with increasing radius moving from E/SO colors in the center to mid- 
type galaxy colors at the limit of the data, r ~ 2 /r'Mpc. The green points 
represent the uncorrected V — i' profile of the green sample. 

for the tendency towards more luminous early-type stars that 
are responsible for the bluer galaxy colors at large radius and 
which otherwise bias the interpretation of the M/L ratio, as 
described below. 

8. CLUSTER LUMINOSITY FUNCTIONS 

The data allow the luminosity function to be usefully con- 
structed in several independent radial and magnitude bins, and 
hence we can examine the form of the luminosity function of 
cluster members as a function of projected distance from the 
cluster center. For this we combine the ACS and the Subaru 
photometry. The ACS has the advantage of extending two 
magnitudes fainter in the /'-band than the Subaru photometry 
for r < 2'. As can be seen in Figure [15] the ACS i' magni- 
tudes agree with Subaru i' magnitudes for galaxies found and 
matched in both catalogs. 

The background correction to evaluate f c i, must be made in 
each magnitude bin independently, since the relative propor- 
tion of background galaxies increases with apparent magni- 
tude, so that the lower luminosity bins of the green luminosity 



FIG. 15. — Comparison between ACS and Subaru photometry for objects 
in common in the central region covered by both datasets, r < 2'. There is 
very good agreement between magnitudes of the two datasets which have 
independent zero-points and of course independeent photometry. 

function are expected to contain a greater fraction of back- 
ground galaxies and hence should have a relatively higher 
value of gj.' . This trend is apparent in Figure [16] (left pan- 
els), where we plot the recovered mean tangential distortion 
(here the average is over a magnitude bin) for each of the four 
radial bins, as a function of absolute magnitude. A clear trend 
is found at all radii towards higher levels of gj at fainter lu- 
minosities. Note that the mean level of the background distor- 
tion (black solid line) drops with increasing radius so that the 

proportion gr(M)/g T is generally an increasing function of 
radius and a decreasing function of luminosity. To correct for 
this we simply apply equation ( [181 ) to each magnitude bin: 

*rfCMt) = • [ 1 - {gf{M k )) I (gf{r)) ] (22) 

(Note that the background signal is averaged over the whole 
range of magnitudes at that radius.) 

We then construct the luminosity function for four indepen- 
dent radjal_binSj_asshown in Figure [16] (middle panel) and 
fit a lSchechterl (jl976) function to each (dashed lines). It can 
be seen that there is no obvious tendency for the shape of 
the luminosity function to change with radius. The faint-end 
slope of a Schechter function fit is a = -1 .05 ± 0.07 in the i'- 
band. This constancy with radius has been argued with some- 
what l ess significance in other well studied massive clusters 
(e.g.. lPracy et alJl2005h . based on similar deep 2-color imag- 
ing, where the limiting radius is more restricted. We also con- 
struct a composite luminosity function for the whole cluster 
(FigfTTIi. for r < 10', which shows clearly the effect of our 
"g-weighted" background correction, without which the faint- 
end slope would be considerably steeper, a ~ 1.4. 

Our approach is of course essentially free of uncertainties 
in the subtraction of background galaxies by its nature. While 
qualitative similarity between the results of the various studies 
is clear, agreement in detail is not necessarily expected, given 
the likely dispersion in the strength of this effect between 
clusters. Also, the question of background contamination is 
always an issue in the standard approach due to the inher- 
ent fluctuations in the surface density of background galaxies, 
and the need to establish the background counts at a suffi- 
ciently large radius from the cluster to avoid self-subtraction 
of the cluster at t he boundarie s of the d ata, a subject explore d 
in depth, e.g. lAdamietalJ (1200 0): Paolillo et aij (1200 ll): 
lAndreon et al.1 i]2005l) : lHansen et all i]2005l) : iPopesso et all 
(120051) . 
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FIG. 16. — Center: Luminosity functions are shown for four independent radial bins, indicating little trend with radius. The left-hand panel is the degree 
of tangential distortion gr as a function of magnitude used in the derivation of the corresponding luminosity function. Notice that in general gj increases with 
decreasing luminosity because the level of background increases at fainter magnitudes. The right-hand panel shows the lcr, 2<r, and 3<r contours for the Schechter 
function parameters M* and a, for each of the corresponding luminosity functions. 

difference between integrating up to a limiting magnitude of 
My < -14 and extrapolating up to M,< < -10 is only about 
0.1%. 

The lack of any obvious upturn in the cluster lumi- 
nosity function to very faint luminosities, M,< < -12, in 
the cluster core, is in agreement with several other stud- 
ies based on deep photometry of large cluster samples, 
e.g., th e comp osite cluster luminosit y function der ived by 
n (199%: iGarilli et alj (fl999h: iPaolillo et alj (l200l) : 
(I2005h : iPracy et al] (l2005h . where wide field 
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FIG. 17. — Composite luminosity function of the cluster (black squares), 
for r < 10', with the fit to a Schechter function displayed in the lower right 
panel. The green circles show the luminosity distribution without the correc- 
tion made for the background dilution. At bright magnitudes, there is little 
difference between the uncorrected and the corrected points, however, at the 
faint end the uncorrected distribution rises, which if uncorrected would over- 
estimate the faint end slope a ~ 1.4 (top right panel), showing clearly the 
magnitude of the background correction, which when accounted for by our 
method results a flat faint end slope, a ~ 1 . 

We also integrate our luminosity functions as a consistency 
check of the luminosity density profiles derived earlier. This 
is done by calculating $ t /(M) in the same radial bins as our 
luminosity profile above, and summing over the magnitude 
bins: 



0.MM,, @ -M k ) 



(23) 



The results shown above in FigureQ~2](blue points) agree very 
well with those of L c \ described in the previous section. 

Note that in constructing these luminosity density profiles 
we have implicitly assumed that the luminosity function is in- 
tegrated over fully. Fortunately our data is complete to a suffi- 
cient depth (j! < 26.5) so that the contribution of the integrated 
luminosity density from undetected objects is very small, as 
evaluated when we examine the luminosity functions. The 
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imaging is employed for several Abell clusters and overdensi- 
ties identified in the SDSS survey, and where careful attention 
can be paid to the background counts and their uncertainty. 
The study of Pracv et al. ( 2005) is the most similar to our own, 
containing three rich and fairly distant Abell clusters, and here 
the LF's show no obvious upturn to M < -13 with a generally 
flat Schechter function slope in the range -1.1 to -1.25. 

For the well studied Coma cluster, a steepe r slope has been 
claimed, a ~ —1.4, bv lBernstein et al.l d 19951). though subse- 
quent faint spectroscopy bv lAdami et al.l d2000l) . has revealed 
the presence of a background cluster at z — 0.5, which when 
corrected for leads to a fiat faint-end slope. In contrast an up- 
turn is clai med for a composite sample of 25 SDSS selected 
clusters by Pop esso et al.l (120051) . though an individual exam- 
ination shows considerable variation, with only a minority of 
~ 6 displaying a distinct upturn which varies in amplitude, so 
that one may wonder about the role of anomalous background 
count fluctuations in these cases. 

9. M/L PROFILES 

We may now go on to examine the mass-to-light ratio using 
the mass profile previously determined based on the distor- 
tion and m agnification profi l e of th e red galaxy sample, as 
derived by Broadhu rst et al.1 (l2005bl) . and also based on the 
central s trong lensing informati on derived from 106 multiple- 
i mages (|Broadhurst et al.ll2005al) . The mass profile derived 
in Broadhurst et al. (2005b) was found to be somewhat more 
pronounced than a simple NFW form, with the observed gra- 
dient increasing monotonically with radius from the cluster 
center. Dividing the mass profile by our newly derived light 
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FIG. 1 8. — Mass-to-ligh t ratio v s. radius. Mass profile is taken from a fit 
made by Broadhurst et al. 1 2005b) using the same data analyzed here (blue 
points). We also plot the total mass-to-stellar mass ratio M/M* (red points), 
accounting for the cluster color profile. 
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FIG. 19. — Cumulative Mass-to-light ratio vs. radius. The profile declines 
markedly inside r < 100 /i~'kpc, otherwise it remains fairly constant with 
radius. 

profile we obtain a profile of the mass-to-light ratio for A 1689 
as shown in Figure [T8l 

We find that the M/L ratio peaks at intermediate radius 
around 2 ' (r~ 100 /r'Kpc) and th en falls off linearly to larger 
radius. (Broadhurst et al. 2005a|) have previously identified 
the drop in M/L at small radius as due principally to the tight 
central clump of very luminous cluster members, noting that 
this clump may have resulted from the effect of dynamical 
friction. 

Note that the peak value is rather large, equivalent to 
M/L B ~ 4OO/z(M/L) in the restframe B-band often used as 
a reference, but the mass of A 1689 is at the extreme end of 
the clu ster population, M~2x 10 I5 /iM Q (Broadhur st et al.l 
2005b|), and given the general tendency of M/L to increase 
with increasing mass, from galaxies through groups to clus- 
ters, we may not be surprised to find the peak to some- 
what exceed the typical range for clusters, 150- 300H.M / Lb q 
(Carlb erg et al .1 12001:!). The general profile o f M/L is similar 
in form to that derived for CL 0152-1357 by lJeeetal.1 (120051) 
based on a careful weak lensing analysis of recent deep ACS 
images. 

We have also constructed a profile of the total mass to stel- 
lar mass M/M* ratio. This is arguably a more physically use- 
ful indicator of the relationship between dark and luminous 



matter compared to the ratio of M/L, because the starlight 
can be strongly influenced by the presence of relatively small 
numbers of luminous hot stars. To calculate M/M* we make 
use of the color profile derived in § [7] and an empirical rela- 
tionship between color and the ratio M*/L for stellar popu- 
lat ions establ i shed f or local galaxies in the the SDSS survey 
by iBell et al.l (|2003). The slope of the projected stellar mass 
profile, dlogM^/dlogr = —1.15 ± 0.13, derived this way is 
slightly steeper than the luminosity profile, as expected. The 
observed relation we derive this way is somewhat flatter than 
for M/L and the mean contribution by mass for stars is about 
1.25% for this cluster and similar to a mean value of ~ 2% 
derived from a carefully selected sample of local clusters by 
iBiviano & Saluccil d2006l) . 

10. CLUSTER MASS PROFILE 

We use the combined distortion information obtained from 
the ACS and Subaru imaging, as described above (Fig. [6} 
and compare with models for the mass distribution. We have 
improved on our earlier distortion measurements made with 
Subaru, with the addition of the background blue galaxy pop- 
ulation defined here, so that the significance of the distor- 
tion measurements is somewhat greater than our earlier work 
which was based only on the red sample (Broadhur st et al.1 
I2005bl) . In addition, we have extended the distortion mea- 
surements to the central region using the HST/ACS informa- 
tion, as described in § |4] where we have clearly identified a 
maximum and a minimum value of gj, which accurately cor- 
respond to the tangential and radial critical curves (Figj20j>, 
independently derived from the the many giant tangential 
and rad ial arcs observed for this cluster (see iBroadhurst et all 
l2005ah . 

Here we test the univers al parameterization of CDM-based 
mass profiles advocated bv lNavarro et al.l dl997l) . This model 
profile is weighted over the differing results from sets of 
haloes identified in N-body simulations. A cluster profile is 
summed over all the mass contained within the main halo, in- 
cluding the galactic haloes. Hence, we compare the integrated 
mass profile we deduced directly with the NFW predictions 
without having to invent a prescription to remove the cluster 
galaxies. 

NFW have shown that massive CDM haloes are predicted 
to be less concentrated with increasing halo mass, a trend 
identified with collapse redshift, which is generally higher for 
smaller haloes following from the steep evolution of the cos- 
mological density of matter. The most massive bound struc- 
tures form later in hierarchical models and therefore clusters 
are anticipated to have a relatively low concentration, quanti- 
fied by the ratio C,,„- = r,,,-™/ jr s . In the context of this model, 
the predicted form of CDM dominated are predicted to fol- 
low a density profile lacking a core, but with a much shal- 
lower central profile (r < 100 h~ l kpc) than a purely isother- 
mal body. 

The fit to an NFW profile is made keeping r s , and p s , the 
characteristic radius and the corresponding density, as free pa- 
rameters. These can be adjusted to normalize the model to the 
observed maximum in the distortion profile at the tangential 
critical radius of ~ 45". The combination of these parameters 
then fixes the degree of concentration, and the corresponding 
lensing distortion profile can then be calculated. 

Integrating the mass along a column, z, where r 2 = 
(£ r r s ) 2 +z 2 gives: 

«<0-A*Q/ (24) 
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Using this mass, a bend-angle of a = 4 ™q^P ^* is produced 
at position 9 = £,r s /d[. The mean interior mass within some 
radius r x = xc x r s can be obtained by integration of the NFW 
profile giving: 



M« r)=—r> 



4tt 3 a„(o)A c (l+z) 3 



x 3 n m (zi) 



\og(l+xc x )--^ 
logd+c,)-^ 



(25) 



where the above integral is carried out to the virial radius. 
Here when we make this comparison we do not attempt to fit 
the distortions in the region of the radial and tangential critical 
curves because the measurements must underestimate the am- 
plitude of the model predictions near these curves due to the 
finite size of the background galaxies, so that the model max- 
imum and minimum, gf = ± 1 corresponding to the tangential 
and radial critical curves cannot be reached by the data (see 
Fig. [6]l near these critical curves. Finite area sources are on 
average not as magnified or distorted as an ideal point source 
due to the gradient of the lensing magnification over the sur- 
face of the source, so that for a source lying over a lens caustic 
only an infinitesimally small area is infinitely magnified. In 
principle simulat ions based on realistic galaxy samples like 
those modeled by lBouwens et al.l (Q998) could be used to cor- 
rect for this effect but this will await further work. 

We can instead utilize here the observed location of these 
critical curves since a clear maximum and minimum is ob- 
served in the inner distortion data (Fig. [6]l and has also 
been independently deter mined from the multiple-image data 
dBroadhurst et al.ll200 5a) for this cluster. In addition to the 
location of the critical curves we can also clearly identify the 
radius where gj = 0, lying in between these critical curves. At 
this radius the radial and tangential magnifications are equal 
and hence images are unchanged in shape (though in gen- 
eral highly magnified), so the observed value of gj will pass 
through zero at a radius in between the critical curves. This 
radius corresponds to the contour where the projected surface 
density is e qual to the critical surface density, £ = £ fn -, (e.g., 
Kaiser l 1995b . and hence is smaller for more concentrated pro- 
files. In Figure l20l we plot these two radii as a function of 
model concentration, where the models are all normalized to 
a tangential critical radius o f 47" to match the mean critical 
radius derived from the data (Broadhurst et al. 20 05al) . 

To normalize the models we choose to reproduce the ob- 
served Einstein radius of 47" and compare the predicted lo- 
cation of the radial critical curve and the k = 1 curve. Fig- 
ure [20] shows that both these radii decrease slowly as the 
concentration parameter is increased. We have marked the 
observed values of these radii as determined by two inde- 
pendent observational means. We can use the statistical 
distortion measurements as described in § [3] and the same 
va lues derived from the multiple image analysis presented 
in Broadhurst et al. (2005a). These differing estimates are 
closely consistent with each other, and by comparison with 
the model curves bracket intermediate values of concentra- 
tion in the range 5 < C l7 > < 15 (Fig. [20l , a range consistent 
with the results from a detailed f it to the inner profile mea- 
sured in Broadhurst et al.l (12005 al) . This is found to be very 
simila r to the independe n tly derived central p r ofile of A 1689 
fromlDiego et al.l d2005l) : IZekser et al.1 d2006b : lHalkola et all 
(120061) . 

Outside the tangential critical radius, for r > 1.5, the dis- 
tortion measurements are small enough not to suffer from any 
significant underestimation due to the finite source sizes and 
we may compare the observed distortion profile, gr(r), out to 



28 

26 

24 

a! 22 

n 

o 

20 



Radial critical curve 
k=1 curve 





FIG. 20. — The curves shows how the radius of the radial critical curve, 
shown in blue, varies with the concentration parameter of an NFW profile, 
where the model is normalized to generate the observed Einstein radius of 
47" for all values of the concentration parameter, C,,,,-. The radial critical 
curve shrinks as the mass profile becomes steeper. This is also the case for 
the radius at which the distortion gr = 0, corresponding to the radius where 
the surface density of matter is equal to the critical surface density (re = 1), 
here the degree of tangential and radial distortion is always equal, indepen- 
dent of the form of the mass profile. We have marked the observed values 
of the radial critical curves and the radius where the distortion is seen to be 
zero as measured independently in two ways using the statistical distortion 
measurements (dotted lines) as described in § f4] and the same va l ues de- 
rived from the multiple image analysis presented in Broadhurst et al. (2005a) 
(dashed lines). These differing measurements are consistent with each other 
and by comparison with the model curves bracket intermediate values of con- 
centration in the range 5 < CW < 15 for the inner strongly lensed region. 



the limit of our data, r < 15". We find this distortion pro- 
file is reasonably well fitted by an NFW profile particularly at 
large radius, 4.0" < r < 15", but with a relatively large con- 
centration, C V ir = 27.2^5 y, as shown in Figure [2]] Note that 
we have used here a linear radial binning when measuring 
the concentration parameter and therefore the result here is 
more weighted by large radius signal than for the analysis of 
Broadhurs t et al.l (1200 5b) where we used logarithmic binning, 
yielding a smaller value of C n >. This difference in the derived 
value of Ci„. is not becuase of any revision in our estimates 
of the distortion, in fact both analyses yield very consistent 
distortion profiles at large radius, but rather that the form of 
the NFW profile is not consistent with our data over the full 
radial range - the best fitting NFW model is either too shallow 
at large radius or too steep at small radius depending where 
one prefers to fit the data. 

We also plot lower concentration profiles, including C t ,,> = 
14, which was found previously by Broadhurst et al. (2005a), 
to fit best the overall lensing derived mass profile from com- 
bining the mass profile derived from the multiply lensed im- 
ages in the central region, r < 2', with the mass distribu- 
tion derived from weak lensing distortion and magnification 
measurements from the red background galaxy sampl e. This 
model fit, as pointed out by iBroadhurst et alj d2005bl) . is not 
as pronounced as the observed surface mass profile, being 
too shallow at larg er radius and too steep at small radius 
(see figures 1&3 of Broad hurst et all l2005bl) . Here we see 
more clearly that this fit with C Vi y = 14 increasingly overpre- 
dicts the observed distortion profile with radius. We also plot 
C V ir = 8.2 which best fit s the central strongly lensed region 
(IBroadhurst et al. 2005a) r < 2', derived from 106 multiply 
lensed images. Again this fit overpredicts the gj profile in 
the we ak lensing regime, as pointed out in IBroadhurst et alj 
(I2005bl) . 
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FIG. 21 . — NFW models are compared with the measured values of gj for 
the red+blue background sample described in §|2] The models are normalized 
to match the observed Einstein radius of 47" . A relatively high concentration 
is preferred with the best fit corresponding to CV = 27+j i (black curve), for 
the bulk of data, r > 1.5', excluding the strong region where the measured 
distortions required a significant correction for the finite source sizes. The 
anticipated low concentration of C ~ 5 (pink curve) is obviously excluded by 
the data. For reference we also overplot the purely isothermal profile normal- 
ized to the observed Einstein radius (red dotted curve). It also overpredicts 
the data. 

We clearly exclude the low concentration profile generally 
predicted by CDM based models of structure formation. A 
value of C v i r ~ 5 is generally anticipated for massive clusters, 
although the s catter in concentratio n at a given mass is con- 
siderable (e.g., Bull ock et al.l l2001). Figure |2TI shows clearly 
how this profile is much too shallow to generate the relatively 
steeply declining observed distortion profile. The triaxial- 
ity of realistic haloes means that projection effects will bias 
somew hat the derived dist orti on profile, as examine d care- 
fully by lOgurietail (120051) and lHennawi etal.1 (120051) . show- 
ing that the level of such bias effect is expected to enhance 
the derived concentration by approximately ~ 20% on aver- 
age. Whilst A1689 is clearly an anomalous cluster in terms 
of the size of the Einstein radius, the cluster is very round in 
terms of the projected X-ray emission, with only minimal sub- 
structure observed in the optical near the center. Hence we are 
left with a clearly unresolved problem, that the observed con- 
centration would seem to far exceed any reasonable estimate. 
Other independent work on the combined profile from strong 
and weak lensin g measurements fo r the clusters C1002 4+17 
(iKneib et alJ2003l) and MS2 137-23 dGavazzi et al.ll2003h also 
point to surprisingly high concentrations, and it is therefore 
important to extend this type of detailed work to other clus- 
ters to test the generality of the profile derived here. 

For reference we also plot the distortion profile for a singu- 
lar isothermal body in Figure [2T| which is simply expressed 
as 

St = ^—t, (26) 

and normalized to the observed Einstein radius, 6>£ = 47". 
This model also overpredicts the data at large radius, indi- 
cating the outer mass profile is steeper than 1 /8 in projection. 

1 1 . DISCUSSION AND CONCLUSIONS 

We have explored a new approach to deriving the luminous 
properties of cluster galaxies by utilizing lensing distortion 
measurements, based on the dilution of the lensing distor- 
tion signal by unlensed cluster members which we assume 



are randomly oriented. We have tested this assumption for a 
restricted sample of bright cluster galaxies which project be- 
yond the faint galaxy background population so the level of 
background contamination is negligible, confirming that the 
cluster galaxies are randomly oriented with a negligible net 
tangential distortion for the purposes of our work. 

This dilution approach is applied to A 1689 to derive the 
radial light profile of the cluster, a color profile and radial lu- 
minosity functions. The light profile is found to be smoothly 
declining and fitted with a power-law slope d log(L) / d log(r) = 
-1.12 ±0.06. We also see a mild color gradient corresponding 
to a change in the cluster population from early- to mid-type 
galaxies in moving from the center out to the limit of our data 
at 2 /r'Mpc. Unlike the light profile the gradient of mass 
profile is continuously steepening, such that the ratio of M/L 
peaks at intermediate radius. We find that the cluster lumi- 
nosity function has a flat faint-end slope of a = -1 .05 ± 0.07, 
nearly independent of radius and with no faint upturn to 
M r < -12. 

A major advantage of our approach is that we do not need to 
define far-field counts for subtracting a background, as in the 
usual method where there is a limitation imposed by the clus- 
tering of the background population that limits the radius to 
which a reliable subtraction of the background can be made. 

We have also established that the bluest galaxies in the field 
of A 1689 lie predominantly in the background, as their radial 
distortion profile follows closely the red galaxies, but with an 
offset indicating the blue population lies at a greater mean dis- 
tance than the red background galaxies, and consistent with 
the estimated mean redshifts of these two populations. With a 
larger sample of clusters, this purely geometric effect can po- 
tentially be put to use to provide a simple model-independent 
measure of the cosmological curvature. 

The mass profile of A1689 was reexamined using our com- 
bined background sample of red and blue galaxies. The dis- 
tortion profile derived from this sample is consistent with our 
earlier work, but somewhat more statistically significant, so 
we have examined the mass profile more carefully out to a 
larger radius. We have found that the distortion profile is 
steeper than predicted for CDM haloes appropriate for cluster 
sized masses C Vi > ~ 5. This discrepancy is particularly clear 
at large radius r > 2', where an acceptable fit is found to an 
NFW profile but with a concentration C n > = 27^y. This find- 
ing is consistent with our earlier work which showed that al- 
though an overall best fit profile of C V! > ~ 14 to the joint strong 
and we ak lensing based data presented in iBroadhurst et alj 
(2005b), the curvature of the data is more pronounced than an 
NFW profile being shallower in the inner region and steeper 
at larger radius, so that the derived value of the concentration 
increases with radius depending on the radial limits being ex- 
amined. 

This result is surprising and may require a significant de- 
parture from the standard CDM model, either in terms of the 
mass content, or the epoch at which the bulk of the cluster was 
assembled. For example, one possibility to achieve earlier 
formation of clusters is to allow deviation from Gaussianity 
of the primordial dens ity fluctuation field , as has been consid- 
ered recently by, e.g.. lSadeh et alj ( 120061) . A1689 is amongst 
the most massive known clusters, and projection effects may 
play a role in boosting somewhat the lensing signal along the 
line of sight. We therefore aim to test the generality of this 
result with a careful study of a statistical sample of clusters. 

Upcoming spatially resolved SZ measurements will add a 
significant new ability to determine cluster mass profiles over 
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a large range of radius, and allow for improved consistency 
checks between the various independent means of estimating 
masses. The combination of X-ray, lensing and SZ measure- 
ments will soon lead to far greater accuracy in understanding 
the nature of cluster mass profiles. 

We plan an improvement to the weak lensing work with 
deeper multi-color imaging from Subaru for measuring reli- 
able photometric redshifts for a sizable fraction of the back- 
ground population. This added dimension of depth will en- 
hance the weak lensing signal and reduce the systematic prob- 
lems of cluster and foreground contamination of the lensing 
signal. We also aim to extend this work to well studied clus- 
ters at lower redshift with archived Subaru imaging and de- 
tailed X-ray and upcoming SZ observations as in principle, 
the lensing signal should be equally strong for lower redshift 
clusters, given the maximal ratio of lens to source distances, 



di s /d s ~ 1, for faint background sources. 
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APPENDIX 

APPENDIX A: CLUSTER GALAXY FRACTION FROM THE LENSING DILUTION EFFECT 

Let us derive eq. ( [T8l l. For simplicity, here we assume the weak lensing limit so that the reduced shear is approximated 
by the gravitational shea r, g sa 7. It is useful to factorize the lensing signal with the geometry-dependent factor such that 
dSeitz & Schneider! 19971) : 



gT (r) = w{z)g T , o{r) 



(Al) 



where gT,oo( r ) denotes the tangential shear calculated for hypothetical sources at an infinite redshift, and w(z) is the lensing 
strength of a source at z relative to a source at z — > 00, w(z) = D(z)/D(z — + 00); D(z) = d[ s /d s as introduced in §5. The relative 
lensing strength vanishes for cluster and foreground galaxies, that is, w(z) = for z < Zi- 
As the tangential shear is obtained by averaging over an annular region, it can be formally written in the following form: 



(grW) =gT,ao(r) 



J d 2 x dz dn / dz w(z) 



= gT,oo(r) 



/ 00 dz dN /dz w(z) 



N t , 



(A2) 



J d 2 x dz dn / dz 

where dn/dz is the surface number density distribution of galaxies per unit redshift interval per steradian, dN /dz = J d 2 x dn/dz 
is the mean redshift distribution of galaxies in the annulus, N lol = J^°dz dN /dz is the total number of galaxies in the annulus, 
Nb g = f^°dz dN /dz is the number of background galaxies in the annulus, (w) z>Zl = J^°dzdN/dz w(z)/ J^dzdN/dz is the mean 
lensing strength without including the dilution effect; here we have assumed that the lensing properties are constant over the 
annulus where we take the ensemble averaging. Note that the factor Nb g /N ror accounts for the dilution effect on the lensing signal 
strength due to contamination by foreground and cluster-member galaxies. In general, there is a contribution from foreground 
galaxies to the total number of galaxies N to t- However, for the case of A1689 at a low redshift of zi = 0.183, this contribution is 
negligible. That is, iV tot i=s N c i +Nb g with N c i being the number of cluster galaxies in the annulus. For a background galaxy sample, 
Ntot =N bg . 

Since we are to compare galaxy samples with different redshift distributions, we need to account for different values of the 
mean lensing strength, (w) z>Zl . As explained above, our green sample (denoted with G) comprises both cluster and background 
galaxies. Hence, according to eq. (1A21 i. the expectation value for the mean tangential shear estimate is 
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As for our background sample (denoted with B), including the red and blue samples, this is 

(g ( T\r))=gT,oo(r)(w ( %.> Zr 
By taking the ratio of the two tangential shear estimates, we obtain the following expression: 



(gf\r))/{gT(r)) 



N 



(G) 



bg 



N-i+Nl 



-(G) 



bg 

Alternatively, we have the expression for the cluster galaxy fraction as 



(By 



Z>Zl 



Mr) i 



N d 



(gf'ir)) (W 



(B)\ 



/z>zi 



N c ,+N 



(C) 
bg 



(gf(r)) (w(«> 



z>z. 



(gf\r)) (0%, 
(#>(r)> (D^) Z>ZI 



(A3) 



(A4) 



(A5) 



(A6) 



This is the desired formula for the cluster galaxy fraction from the weak lensing dilution effect. In order to take into account 
different populations of background galaxies in the two samples, one needs to estimate the correction factor, (D (B) ) z>Zl / (D (G) ) ,>,, . 
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APPENDIX B: NON-LINEAR EFFECT IN THE REDUCED SHEAR ESTIMATE 

In Appendix A, we assume that the observable reduced shear is linearly proportional to the lensing strength factor, w(z). 
However, the reduced sear, defined as g = 7/(1 — k), is non-linear in k, so that the averaging operator with respect to the redshift 
generally acts non-linearly on the redshift-dependent components in g. 

To see this effect, we expand the reduced shear with respect to the convergence k as 

00 

g= 7(1 -K)" 1 = WJooil-WKoo)' 1 = WJoo ^(vVKoo)' (Bl) 

k=0 

where and 700 are the lensing convergence and the gravitational shear, respectively, calculated for a hypothetical source at an 
infinite redshift. Hence, the reduced shear averaged over the source redshift distribution is expressed as 



00 

(S)=7oo5>* + Voo- (B2) 

k=0 



In the weak lensing limit where Koc, \j\oo <^ 1> then (g) w (w)7oo- Thus , the mean reduced shear is simply proportional to the 
mean lensing strength, (w) . The next higher-order approximation for eq. (IB2b is given by 



(g) * 7oo «w) + (w 2 )^) « ; {W ) 1 7, I , y (B3) 

l-/too(W z )/(w) 

ISeitz & Schneider! (Q997) found that eq. ( |B3t yields an excellent approximation in the mildly non-linear regime of < 0.6. 
Defining f w = (w 2 ) / (w) 2 , we have the following expression for the mean reduced shear valid in the mildly non-linear regime: 

(g) * T^Tua W 

with (k) = (w)k 00 and (7) = (w)7oo dSeitz & Schneiderl! 1997b . For lensing clusters located at low redshifts of zi < 0.2, (w 2 ) ~ 
(w) 2 or /,„ « 1, so that (g) « ( 7 )/(l - <«)). 

The ratio of tangential shear estimates using two different populations B and G of background galaxies, in the mildly non-linear 
regime, is given as 

The lowest-order correction term is proportional to (/, ( t B) -/iP(w (G) )) Kqq, which is much smaller than unity for the galaxy 
samples of our concern in the mildly non-linear regime. In conclusion, it is therefore a fair approximation to use eq. ( fTHT l for 
measuring the cluster galaxy fraction via the dilution effect. 
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